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Abstract 

The QCD vacuum condensates and various vacuum susceptibilities are all important parameters which char¬ 
acterize the nonperturbative properties of the QCD vacuum. In the QCD sum rules external field formula, 
various QCD vacuum susceptibilities play important roles in determining the properties of hadrons. In this 
paper, we review the recent progress in studies of vacuum susceptibilities together with their applications to 
the chiral phase transition of QCD. The results of the tensor, the vector, the axial-vector, the scalar, and the 
pseudo-scalar vacuum susceptibilities are shown in detail in the framework of Dyson-Schwinger equations. 
Keywords: vacuum susceptibility, Dyson-Schwinger equations, chiral phase transition 


* Corresponding author 

Email addresses: phycuiOnju.edu.cn (Zhu-Fang Cui), zonghs@nju.edu.cn (Hong-Shi Zong) 


Preprint submitted to Annals of Physics 


May 5, 2015 





Contents 


1 Introdnction 3 

2 Derivation of the vector vacuum susceptibility 6 

3 Calculations of various vacuum susceptibilities 9 

3.1 Calculations of quark propagator and quark-meson vertex. 9 

3.2 Various vacuum susceptibilities. 14 

3.2.1 The tensor vacuum susceptibility. 14 

3.2.2 The vector and axial-vector vacuum susceptibilities. 16 

3.2.3 The scalar vacuum susceptibility. 21 

3.2.4 The pseudo-scalar vacuum susceptibility. 27 

4 Susceptibilities and the chiral phase transition of QCD 30 

5 Summary 35 


2 








1. Introduction 


Nowadays, the commonly accepted theory that describes the strong interaction is Quantum Chromody¬ 
namics (QCD). Due to the asymptotic freedom nature of QCD, in the high energy and high momentum trans¬ 
fer region reliable calculations can be done using the perturbative theory. Numerous comparisons between 
the results from theoretical calculations and experimental measurements have shown that QCD is correct in 
the high energy held. However, in the low energy and low momentum transfer region the coupling constant 
as of QCD becomes large and running, consequently things have to be treated non-perturbatively; and 
thanks to the non-perturbative nature, the QCD vacuum is not trivial. People believe that the non-trivial 
properties of QCD vacuum are closely linked to spontaneous chiral symmetry breaking and conhnement, 
which are the two main features of low-energy hadron physics. In order to characterize the QCD vacuum, 
people often introduce various vacuum condensates (such as the two-quark condensate, the gluon conden¬ 
sate, the mixed quark gluon condensates, and four-quark condensate, etc.) and many vacuum susceptibilities 
(such as the tensor vacuum susceptibility, the vector and the axial-vector vacuum susceptibilities, the scalar 
and the pseudo-scalar vacuum susceptibilities) phenomenologically [1, 2, 3, 4, 5, 6]. In this paper, we review 
the recent progress in studies of vacuum susceptibilities together with their applications to the chiral phase 
transition of QCD. The results of various vacuum susceptibilities are shown in detail in the framework of 
Dyson-Schwinger equations (DSEs). 

Metaphorically speaking, we can regard the QCD system as a black box, one way to study the QCD 
vacuum is to cause a perturbation by adding an external field to the system, such that we can learn something 
about the QCD vacuum indirectly by studying its responses to the external held.^ In the theory of QCD, 
the quark propagator is the simplest Green function (two-point Green function), and its linear responses 
to the external helds reflect the properties of QCD vacuum directly. QCD vacuum susceptibilities are 
just the parameters that are associated with the linear responses of the quark propagator to the external 
fields, which can then characterize the non-perturbative properties of QCD vacuum. In the QCD sum 
rules external field formula, QCD vacuum susceptibilities play important roles in determining the hadron 
properties [2, 3, 5, 6]. For example, the strong and parity-violating pion-nucleon coupling depends crucially 
upon the vacuum susceptibilities of pion [7, 8] , and the tensor vacuum susceptibility is closely related to the 
tensor charge of nucleon [9, 10, 11, 12, 13, 14, 15, 16, 17]. The nonlinear susceptibilities are also found to 
be correlated with the cumulant of baryon-number fluctuations in experiments [18, 19, 20, 21, 22]. Due to 
the importance of QCD vacuum susceptibilities in low energy hadron physics field, many different methods 
and models to study various QCD vacuum susceptibilities have been utilized over the years. Take the scalar 
vacuum susceptibility as an example, there are Lattice QCD calculations [23, 24, 25, 26], calculations using 
the DSEs method [27, 28, 29], the multi-flavor Schwinger model [30, 31], the linear sigma model [32, 33], 
the (Polyakov loop extended) Nambu-Jona-Lasinio ((P)NJL) model [34, 35, 36], the nonlocal chiral quark 


^In the experimental side, the physicists often measure the linear responses that are proportional to the perturbations, such 
as many susceptibilities, conductivity, etc. Therefore, in order to compare with experimental results the calculations of various 
linear responses are very important in the theoretical side. 
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model [37], and the chiral perturbation theory [38, 39], etc. However, owing to the lack of a model independent 
expression for calculating the vacuum susceptibilities, different models or approximations often give different 
results, some of which may even be distinct. In this case, in order to provide a good starting point for 
studies of vacuum susceptibilities, the authors of Ref. [40] proposed a model independent equation for 
calculating the vector vacuum susceptibility firstly, using the QCD sum rules external field formula. It 
should be pointed out that, their results are based on the linear response theory to an external field of 
fermions, so that it is universal, and can also be used in other quantum Held theories (e.g., Quantum 
Electrodynamics in (2+1) dimensions, QED 3 ) [41, 42, 43]. Based on this, people have also studied the 
axial-vector and the tensor vacnnm susceptibilities in the framework of DSEs and under the rainbow-ladder 
approximation [44, 45], while the authors of Ref. [46] studied the scalar vacuum susceptibility beyond the 
rainbow-ladder approximation, as well as the pseudo-scalar vacuum susceptibility [47]. Here we want to stress 
that the authors of Ref. [47] derived a model independent result for the pseudo-scalar vacuum susceptibility 
using the isovector-pseudo-scalar vacuum polarization, whereas the authors of Ref. [48] obtained model 
independent results for the vector and axial-vector vacuum susceptibilities using the vector and axial-vector 
Ward-Takahashi identities (WTIs). In Ref. [49], the authors calculated the tensor vacuum susceptibility 
employing the Ball-Chiu (BC) vertex [50, 51], and also studied the dressing effect of the quark-gluon vertex 
on the tensor vacuum susceptibility. Their results show that the tensor vacuum susceptibility obtained in 
the BC vertex approximation is reduced by about 10% compared to its rainbow-ladder approximation value, 
which means that the dressing effect of the quark-gluon vertex is not very large in the calculation of the 
tensor vacuum susceptibility within the DSEs framework. These works promoted our understandings of the 
five kinds of vacuum susceptibilities, and up to now, only the scalar and the tenser vacuum susceptibilities 
are still model-dependent, to choose a more reliable model then becomes very important. 

It is commonly accepted that with increasing temperature and/or quark chemical potential, strongly 
interacting matter will undergo a phase transition from the hadronic phase to the quark-gluon plasma 
(QGP) phase, which is also expected to appear in the ultra-relativistic heavy ion collisions [52], and the 
chiral symmetry which has been broken dynamically is also supposed to restore partially [53]. This is a very 
important as well as a very hot topic. One of the main goals for RHIC (Relativistic Heavy Ion Collider) and 
LHC (Large Hadron Collider) is just to create and then study such a new state of matter. These studies 
also play crucial roles in researches on the evolution of the early universe. The properties of the QCD 
vacuum will change when the temperature and/or quark chemical potential vary, as a result the vacuum 
condensates and vacuum susceptibilities which reflect the nature of the QCD vacuum should also change with 
temperature and/or quark chemical potential. Theoretically, the QCD phase transitions via the calculations 
of temperature and quark chemical potential dependence of the two-quark condensate or QCD vacuum 
susceptibilities can then be studied. In the standard definition, an order parameter is a quantity which is zero 
at one side of a phase transition and nonzero at the other side, so that the two-quark condensate in the chiral 
limit is such a quantity for the chiral phase transition (therefore, in some cases it is also referred to as “chiral 
condensate”) [54]. However, while away from the chiral limit the two-quark condensate is no longer an exact 
order parameter (since there is no exact phase transition), although it can be stretched to encompass first 
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order transitions where the quantity has jumps from one nonzero value to another. Various susceptibilities, 
which are related to the derivatives of the two-quark condensate, are usually used to locate the transition, 
and hence are also important quantities to characterize the properties of the system.^ For example, the scalar 
susceptibility is often used to describe the QCD chiral symmetry restoration at finite temperature and finite 
chemical potential [23, 28, 35, 55, 56, 57, 58, 59, 60], while studies of the temperature and chemical potential 
dependence of the quark number susceptibility can provide useful information to the vicinity of the critical 
end point (CEP), where the first order phase transition meets with the crossover and a second order phase 
transition takes place [61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 46, 79, 80, 81, 82], 
The cases with finite temperature and chemical potential represents a big branch in the field of studies 
on QCD vacuum susceptibilities. Hence it is difficult to give a comprehensive presentation on its theories 
and applications to the chiral phase transition of QCD. In this paper, we mainly discuss the calculations 
of various QCD vacuum susceptibilities at zero temperature and zero chemical potential, and then present 
some results in studies of the chiral phase transition of QCD. We will review the QCD sum rules external 
field formula in a nutshell, then take the constant vector field as an example to give the general derivation 
of the vector vacuum susceptibility, and obtain a model independent expression at last. We will also discuss 
the tensor, the vector, the axial-vector, the scalar, and the pseudo-scalar vacuum susceptibilities in detail 
together with their applications. From the model independent form of the vacuum susceptibilities, we 
see that their calculations are closely linked to the quark propagator (two-point Green function) and the 
corresponding vertex function (three-point Green function). Lattice QCD, which is commonly accepted 
as a first principle calculation of QCD itself, can give good results for the two-point Green functions, but 
for three- and more points Green functions there is still a long way to go. Accordingly, in order to do 
the calculations we have to, at present, turn to some effective models or approaches that are based on 
QCD itself. In this paper we mainly use the framework of DSEs approach, which is a very useful and 
successful non-perturbative method to treat the non-perturbative strong interactions as well as the properties 
of hadrons [83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96]. It must be emphasized here that, DSEs 
is an infinite series of coupled non-perturbative integral equations that the Green functions of continuous 
space-time should satisfy, of which a given n-point Green function is related with the m-point Green function, 
where m > n. In order to solve practical problems, people have to take “appropriate” approximations to 
the DSEs.^ So far, the most commonly used approximation scheme is the rainbow-ladder approximation. 
The rainbow approximation means that in the DSE for quark propagator the full quark-gluon vertex is 
replaced by the bare vertex, while the ladder approximation is that in the Bethe-Salpeter equation (BSE) of 
the quark-meson vertex the quark-anti-quark scattering kernel is replaced by its lowest order perturbative 
contribution, details of which can be found in Refs. [83, 85, 88]. The rainbow-ladder approximation is self 
consistent [97].^ Besides, rainbow-ladder approximation also promises the DSEs to satisfy the axis vector 

^Moreover, the two-quark condensate is actually unmeasurable, while various susceptibilities are reachable in the experiments, 
good approximation should holds the original symmetries of QCD itself as many as possible, and meanwhile, easy to 

calculate. 

^In the global color symmetry model (GCM)[86, 98], if adopting the saddle point approximation, it can be found that the 
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WTI, which ensures that pion is the Goidstone boson, and its mass is zero in the chirai iimit. Through many 


years of research, it has been found that the DSEs under rainbow-iadder approximation have achieved great 


success in describing the pseudo-scaiar and the vector mesons, since in these cases the contributions from 
higher order correction terms is iittie [85, 87]. However, the resuits for the scaiar meson are not satisfactory, 
since in this case the higher order corrections contribute significantiy [84]. Therefore, for the scaiar vertex 


and the scaiar vacuum susceptibiiity the rainbow-iadder is not a good approximation. In other words, the 


resuits obtained in this approximation are not reiiabie, and one must go beyond it [46]. Accordingiy, in 
the foiiowing discussions the tensor vacuum susceptibiiity is within the rainbow-iadder approximation (the 
resuits for the BC vertex can be found in Ref. [49], with a decrease of about 10%), whiie for the scaiar 
vacuum susceptibiiity we try to go beyond it. 

2. Derivation of the vector vacuum susceptibility using the QCD sum rules external field 


formula 


In this section we derive a model independent expression of the vector vacuum susceptibility using the 
QCD sum rules external field formula. To be specific, we will work in Euclidean space throughout this paper. 
The vacuum susceptibilities are closely related to the linear responses of the full quark propagator to the 
external fields, so we add to the Lagrangian a term that the quark current coupled to an external field, 


= qiy)^qiy)Vr{y) 


( 1 ) 


here J^{y) is the quark current, Vr(j/) is the external field, q{y) is the quark field, T is a matrix in the 
direct product space of Dirac, flavor, and color. The linear response term of the full quark propagator to 
the external held can be written as [40], 



( 2 ) 


where (cc') and (a/3) are the color and spinor indices of the quark, jO) is the exact vacuum; Gq’^^{x) and 
Gg’^^{x) denote the non-perturbative and perturbative effects of the full quark propagator G’^p^q{x) under 
an external held, respectively. In the paper we adopt such a concept that, G^p^q{x) is calculated in the 
Nambu vacuum, of which dynamical chiral symmetry breaking (DCSB) and quark conhnement are two of 
the most important features; while t/g’^^(x) is calculated in the Wigner vacuum, where the system holds 
chiral symmetry (in the chiral limit case) and the dressed quarks are not conhned (more discussions can be 
found in Sec. 3.1 of this paper, and Refs. [40, 83, 88]). As we will show in the following part, Gq’^^{x) is 
related to the vacuum susceptibilities and hard to calculate directly, while G^js^qix) and Gq’^^ix) can be 
calculated, and then we can get Gq’^^{x) as well as the vacuum susceptibilities indirectly. 


In the QCD sum rules theory, the vacuum susceptibility is defined as [2, 3, 5, 8], 


^ (6l:g^(x)g-^'(0):l0)^r 

= (rVr)„^ 6,,,x^H{x){0\ : g(0)g(0) : JO), 


( 3 ) 


equation which quark propagator satisfy is just the unrenormalized DSE under rainbow approximation, while the one which 
quark-meson vertex satisfy is just the unrenormalized BSE under ladder approximation. 
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|0) represents the perturbative vacuum, is the vacuum susceptibility, while H{x) is a function that 
characterize the non-local property of the non-local two-quark condensate, and H(0) = 1. From Eq. (3) 
we can see that, in order to calculate we should know {x) first, which means we should know 

Qq’^'^{x) and Q^{x) simultaneously. Therefore, how to solve Oq-^'^{x) and Gq{x) self-consistently becomes 
the key to get a model independent result. Next we will take a constant vector field as an example to 
show how to get a model independent expression for the vector vacuum susceptibility [40] . 

When there is an external field, a term which is the coupling of the vector current and the vector external 
field should be added to the action of the system, namely, AS = f d^x g(x)ijf_ig(x)Vf^(x). In this case, the 
dressed quark propagator in the chiral limit is 

g[V]{x) = J VgVgVA g(x)q(0) exp |-S'[g, q, A] - J d'^x g{x)x^g{x)Vf,{x)^ , (4) 


where 

S[g,g,A] = J , (5) 

and — d^A^ + gf°'^‘^A^^Al. Here we do not show the gauge hxing term, as well as the ghost field 

and its integration measure. Considering the linear response of the quark propagator to the external field, 
we can obtain 


g\y]{x) = J VgVgVA [( 7 (x)g( 0 )] exp {-^[g, g, H]} 


- / VgVgVA 


gix)g{0) / d y g{y)i^,g{y)V^{x) 


exp{-S[g,g, A]} -b ■ 


= (0|rr[g(a::)q(0)]|0) - / d‘^y{0\Tr[g{x)g{0)g{y)'yt,g{y)]\0)Vf,{y) 


( 6 ) 


Tt represents the time-ordering operation on Euclidean time r, G{x) = (6 |Tt- [q(a;)g(0)] |0) = g[V = 0](x) is 
the full quark propagator without an external field (V^ = 0), while 


g'^ix) = (0|T.r[g(a;)g(0)]6)v = - J d‘^y{0\TT[gix)g{0)g{y)xt,g{,yW)'^t^iy) 


(7) 


is the linear response term of the quark propagator to the external field. Hence, now the main task is to 
know f d'^y(0lTrlg(x)g(0)g(y)j/j,g(y)]l0}, we can try to do Taylor expansion of the inverse of dressed quark 
propagator to [40, 44, 45], 


g[v]-^ = g[v]- 


sg-^[v] 


lv„=0 




Vu + ■ ■ ■ = G ^ -b r ^ -b 


v „=0 


and then only keep the first order term of the external field, 


( 8 ) 


g[V] = G-GTf,Vf.G + ---. 


(9) 


Here the vector vertex T^ is defined as 

Tf,{yi,y2-,z) 


Sg[V]{yi,y2) ^ 

6V^{z) 


( 10 ) 
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The derivation of Eq. (9) is general, so is applicable to both full and perturbative quark propagators. In 
configuration space it can be written explicitly as 


C/[V](a;) = G{x) - J d'^ui J d'^U 2 G{x,ui)Tfj^{ui,U2)G{u2,0)Vfj,-\ - 

= G{x)- I dV I d'^U 2 G{x,m) I G{u2,0)V^ 


= G[x)-j 


d^p 


dP--G{p)T^[p,Q)G{p)V^ 


Then the linear response term of the full quark propagator {x) to an external field is, 


g^{x) = (o|r, [g(x)d(o)] |o)v = - 


d'^p 


d’^--G{p)T^{p,0)G{p)V^. 


( 11 ) 


( 12 ) 


and similarly for the perturbative quark propagator (x) 


-•V,PT 


(x) = (0|T. [g(xMO)] |0)z 

= — J d^ui J d‘^U2 G^^{x,ui) 

f d'^P 


(27r) 




( 13 ) 


where G^^{x) = (0|r.r [9(2;)d(0)] |0)i and r^’^(p;0) is the corresponding perturbative vertex. 

Now let’s turn to the calculation of vacuum susceptibilities, the definition of the vector vacuum suscep¬ 
tibility is, 

(6| : g“(0)g^(0) : lO)^ ' V = -^(t • VU5a,X^{b\ : q{Q)q{0) : |0), (14) 

where = (6| : 9a(0)9^(0) ^ |6)v'^ is the local two-quark condensate. It can be seen that the 

two-quark condensate is an explicitly gauge-invariant quantity, and is indeed computed in a gauge-invariant 
way for example in Lattice QCD simulations. That it depends on the renormalization scheme is on the other 
hand obvious.^ If we multiply SabXu°‘ both sides of Eq. (14), then 


(6| : g(0)7.g(0) : \0)^^ ■ V = (-i)x'^V,(6| : g(0)g(0) : |6) 

= [(6|r [g(0)7.g(0)] |6)v - (0|T [g(0)7.g(0)] |0)(;^] • V. (15) 

From Eqs. (12) and (13) we can find that 


[(0|T[g-(0)7.g(0)]|6)v] 

•V = trDcbuGT-VG], 

(16) 

[(0|r[g(0)7.g(0)]|0)(;n 

-V = trocb.G^^r^^-VG^^]. 

(17) 


^Without truncation, DSEs is equivalent to QCD; However, it is impossible to solve infinite coupled equations, to choose a 
better truncation scheme is then not only necessary but also challenging. Concerning the choice of gauge, we will use Landau 
gauge throughout this paper, which has many advantages [83, 88, 95, 96, 99, 100, 101, 102, 103], for example, it is a fixed 
point of the renormalization group; that gauge for which sensitivity to model-dependent differences between Ansatze for the 
fermion-gauge boson vertex are least noticeable; and a covariant gauge, which is readily implemented in simulations of Lattice 
regularized QCD (see, e.g.. Refs. [104, 105, 106] and citations therein and thereto). Importantly, capitalisation on the gauge 
covariance of Schwinger functions obviates any question about the gauge dependence of gauge invariant quantities. 













substituting them into Eq. (15), we can get a QCD sum rules based, also model independent expression of 
the vector vacuum susceptibility as, 

X^(6| : g-(0)g(0) : |0) = ^ [trDcb.GT^G] - (18) 

Here we note that, from the viewpoint of Feynman diagram, trocbuOT^G] and trucbvG^'^T ^'^are 
respectively the values of the full and the perturbative vector vacuum polarizations at zero momentum, so 
the calculation of the vector vacuum susceptibility comes down to the calculation of the vector vacuum 
polarization. Here we want to stress that, as pointed out in Ref. [40], the derivation in this section is general 
and model independent for a constant external field, and so it can also be used to derive other vacuum 
susceptibilities similarly, such as the tensor vacuum susceptibility. However, for the axial-vector vacuum 
susceptibility, we need to take a variable external field, details can be found in Ref. [44]. The signihcance 
of the works in Refs. [40, 44] is that they give theoretically model independent results of various vacuum 
susceptibilities, which provides a good starting point for further model calculations. In the next section we 
will discuss the tensor, the vector, the axial-vector, the scalar, and the pseudo-scalar vacuum susceptibilities 
respectively. 

3. Calculations of various vacuum susceptibilities 

After the model independent expressions of QCD vacuum susceptibilities are obtained, the next step 
is to use these expressions to calculate the values of these vacuum susceptibilities. From Eq. (18) we 
know that, in order to calculate the vector vacuum susceptibility, we must first know the quark propagator 
G{P) (two-point Green function) and the vertex function F^(P, 0)(three-point Green function). This holds 
for the other vacuum susceptibilities. However, as we mentioned earlier, nowadays it is difficult or even 
impossible to know these Green functions from the first principles of QGD directly, and people have to resort 
to various non-perturbative approaches besides the Lattice QGD, such as the QCD sum rules, the chiral 
perturbation theory, the DSEs theory, etc. Generally speaking, each approach has its own advantages and 
disadvantages. Since the early ninety’s of the last century, the DSEs has made great progress, and also 
gained a lot of results which are in good agreement with the experimental measurements and Lattice QCD 
calculations [83, 88, 89, 90, 91, 92, 93, 94, 95]. This approach, therefore, is getting more and more attentions 
in many fields. In the following part of this section, we will discuss the calculations of the quark propagator 
and the quark-meson vertex in the framework of DSEs, and then discuss the calculations of several vacuum 
susceptibilities. 

3.1. Calculations of quark propagator and quark-meson vertex in the framework of the DSEs 

Using the functional path integral techniques, we can get the DSEs for Green functions from the La- 
grangian (details can be found in the textbooks of Quantum Field Theory). Since the derivation of functional 
path integral does not rely on the perturbative theory, DSEs is strict (in the weak coupling limit, DSEs give 
the standard Feynman diagram expansion of the Green functions). The DSE of the renormalized quark 
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propagator is 


G{p) ^ = Z 2 {i'r ■ p + mtm) + Zi q)'y^,^G{q)Tl{q,p), (19) 

here glD^^{p) is the renormalized gluon propagator, A“ is the Gell-Mann matrices with a the color index, 
T'^{q,p) is the renormalized quark-gluon vertex, mhm is the bare mass of current quark and is related to 
A, = f^d‘^q/(27r)^ represents some regularization scheme that keeps translation invariance, where A is 
the mass scale, we should take the limit A —>■ oo at last. A^) and Z 2 (^^, A^) are the renormalization 

constants of the quark-gluon vertex and quark wave function, which depend on both the renormalization 
point p and the renormalization mass scale A. According to the Lorentz structure analysis, the solution of 
the quark DSE Eq. (19) has the following form 

G(p)-i = *7 • pA{p\p^) + B{p\p^). (20) 

Eq.(19) must be solved under certain renormalization condition. Since the asymptotic freedom nature of 
QCD, people usually demand that at a big space-like momentum square there should be 

G{jp)~^ Ip2=^2 =il-p + m(/r), ( 21 ) 

where m{p) is the renormalized current quark mass at p. When chiral symmetry is explicitly broken, the 
renormalized mass and the bare mass satisfy this relation, rn{p) = m(,m(A)/Zm(/r^, A^), here Zm{p\A^) 
is the renormalization constant of mass; while there is no explicit chiral symmetry breaking, Z 2 mbm = 0 , 
which is just the case of chiral limit. Most of the following discussions on the vacuum susceptibilities are for 
the chiral limit case. 

It can be seen clearly from Eq. (19) that, the quark propagator, the gluon propagator, and the quark- 
gluon vertex are related to each other. The gluon propagator and the quark-gluon vertex also satisfy their 
own DSEs, which are coupled to the Green functions with more points. Therefore, in order to obtain a 
closed equation of the quark propagator, we must make some appropriate approximations to the gluon 
propagator and the quark-gluon vertex. On the gluon propagator side, the researches of its DSE have made 
great progress, and the gluon propagator in Landau gauge is for all phenomenological purposes sufficiently 
well-known [83, 88 , 95, 96, 101, 99, 102, 100, 107, 108, 109]. The lack of knowledge refers much more to 
the quark-gluon vertex. It is generally believed that, in studies of the color conhnement, the DGSB, and 
the properties of bound states, the non-perturbative quark-gluon vertex plays an important role [ 110 , 111 , 
112, 113, 114, 115, 116]. One of the most complete investigations of the quark-gluon vertex can be found 
in Ref. [100], and further guides as well as some important recent works and studies beyond the rainbow- 
ladder truncation can be found in Refs. [117, 118, 119] (and references therein). Therefore, when studying the 
dynamical chiral symmetry breaking problems of the quark propagator, people usually choose an appropriate 
gluon propagator model as input. A gluon propagator model often contains some parameters, which can be 
hxed by fitting the low-energy hadron physics experiments (such as the mass and decay constant of pion). 
For the quark-gluon vertex, the case is more complicated. The most simple approximation is just the rainbow 
approximation, namely, use the bare quark-gluon vertex 7 ^ to replace the full one r^(g,p). This is a big 
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approximation, since the momentum dependence of the quark-gluon vertex is completely neglected, and of 
course, the non-perturbative information is also lost completely. In studies of low-energy hadron physics, 
such as the calculation of hadron masses, decay constants, and other physical quantities, it is necessary to 
solve the BSE that the quark-meson vertex satisfy. As mentioned above, the rainbow approximation of the 
quark DSE and the ladder approximation of the quark-meson vertex BSE are consistent with each other, 
and over many years of studies it has been demonstrated that although the rainbow-ladder approximation 
is the most simple one, it can give successful descriptions to some meson properties. The reason why it can 
not describe the scalar meson well has also become known. There have been many attempts that try to 
study the scalar meson beyond the rainbow-ladder approximation. However, this is rather difficult work, 
since if we do some approximations to the quark-gluon vertex of the quark DSE,® we also need to find 
a corresponding quark-meson BSE that is consistent with these approximations, which is really a difficult 
work. As we mentioned earlier, for studies of scalar vacuum susceptibility people should go beyond the 
rainbow-ladder approximation. The authors of Ref. [46] noted that in the calculations of the scalar vacuum 
susceptibility, only the scalar vertex with zero momentum exchange r(p, 0) is used, rather than the full scalar 
vertex T{p,q). Using this feature and the WTI, they give a method that can describe the scalar vacuum 
susceptibility beyond the rainbow-ladder approximation. In the following we will introduce the equations of 
the quark propagator and the quark-meson vertex under rainbow-ladder approximation first, then present a 
method to go beyond it. 

Under the rainbow approximation and taking the chiral limit {m-bm = 0), and substituting Eq. (20) into 
Eq. (19), it is easy to know the equations that A{p^,pL^) and satisfy 


[A{p\n^)-Z2]p^ = - 


glDjjp- 


p-q + 2 


p- (p-q) q-{p-q) 


[p - qf 


B{p\p-^)=A f 

J a 


gsD{p- q)B{q^,p^) 


q^A^{q^,p^)+B^q^,p^)’ 

where we take the Landau gauge for the gluon propagator D“®(p) = (5“® II(p^). After 

inputting the model gluon propagator, Eq. (22) can be solved by numerical iteration. From Eq. (22) we 
can see clearly that it has two distinct solutions: the Nambu-Goldstone solution (or Nambu solution) that 
^ 0 and the Wigner-Weyl solution (or Wigner solution) that B{p^^p?) = 0. In the Nambu- 
Goldstone phase: 1) chiral symmetry breaking occurs dynamically, the originally massless current quarks 
obtain masses by DCSB; 2) quarks are confined, since the quark propagator does not have the Lehmann 
representation. While in the Wigner phase, there is no DGSB, and the quarks are not confined [83, 86, 88]. 
In the Wigner phase, Eq. (22) becomes. 


[A'{p\p‘^)-Z'^\p^ = - 


glDjjp - q) 

q^A'{q\p^) 


p-q + 2 


p- {p-q) q- {p-q) 

{p-q)^ 


(23) 


®For example, by the constraint of the longitudinal WTI and the requirement that the vertex has no kinematic singularity, 
people have developed the Ball-Chiu (BC) and Curtis-Pennington (CP) approximations for fermion-boson vertex in the gauge 
theory [50, 51, 120]. In recent years, people also try to use the constraint of the transverse WTI [121, 122]. 
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here represents the vector part of the self-energy function. Therefore, the chiral limit case of the 

quark propagator in the Wigner phase is 




—17 • p 
A'{p^, pL^)p‘^' 


(24) 


In the ladder approximation, quark-meson vertex (three-point Green function) also satisfies its own BSE. 
For the tensor vertex, its BSE is 

r-(p,0) = 9)7,G(g)r-(g,0)G(g)7c, (25) 

where the superscript m denote “meson”, and then E means quark-gluon vertex while E™ represents quark- 
meson vertex. Similar to the quark propagator, we choose such a the renormalization condition 




(26) 


According to Lorentz structure analysis, the general form of can be written as the following 


rjTi/b^O) = ,p) + a^,^p^p^K 2 {j) ,p)+if,PvA'i{p ,p) 

+7 • PPp.Py^iip^, + iPtJ.P,^A5{p'^, p^). 


(27) 


Now substituting Eq. (27) into Eq. (25) and after some algebra, we get five coupled integral equations of the 
scalar function Ai(p^, /r^)(i = 1, • • •, 5), which have the solution A 3 ^ 4 _ 5 (p^, p?) = 0. So that [45] 

r]rj/(7')0) = a^^Ki{p^, y?) + (j^^p^p^K2{p^, y?), (28) 


The BSE that the vector vertex satisfy is 

E;r(p,0) = Z27^ j\lD^^^p-q)^r,G{q)T^{q,Q)G{q)^^. (29) 

Here one way to solve this is to write down the general Lorentz structure of E™(p, 0) and then do the numerical 
iteration, another more direct way is to use the Ward identity (which is obtained via the Ward-Takahashi 
identity by taking the infrared limit of the photon) 


rir(p,0) = - 


dG{p)- 

dpa 


(30) 


by calculating the partial derivative of the quark propagator to the momentum p, we can know the vector 
vertex directly. Similar treatment can also be used to the axial-vector vertex, for the details please see 
Ref. [44]. 

For the scalar vertex, the rainbow-ladder approximation is not a good one and we have to go beyond it. 
Let us see its BSE first. 


E™(p,0) = l + I ^K{q,p)[G{q)r^{q,0)G{q)], (31) 

here K{q,p) is the quark-anti-quark scattering kernel, which is a four-point Green function that contains 
many non-perturbative information. As mentioned above, when beyond the rainbow-ladder approximation it 
is very difficult to find the corresponding BSE of the vertex. But if take a closer look we can find that only the 
scalar vertex (it is known from the functional path integral theory that taking the derivative of the two-point 
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Green function with respect to the scalar external field will give the scalar vertex) 0) = —|m=o 

is needed to calculate the scalar vacuum susceptibility. The authors of Ref. [46] skillfully used this point. 
By the DSE of quark propagator, Eq. (19), it is easy to know the equation of scalar vertex 

r-(p,o) = i + Ei(p) + S2(p), 

here 

Si(p) = -jj\^D^,{p-q)^^G{q)r^{q,0)G{q)r,{q,p), 

^2{p) = ^ g^D^^{p-q)')^G{q)K{q,p), (32) 

where the summation over color is carried out to be 4/3 and Ky(q^p) = ■ 

Using the vector Ward-Takahashi identity,^ 

ikyT„(q,p) = G{q)~^ - G(p)~^, k = q-p (33) 

we can know another identity that A^(q,p) and r^(p) satisfy, 

ik,A,{q,p)=r^(q,0)-T^{p,0). (34) 

This is an important step towards the consistent solutions of quark propagator and the scalar vertex. 


By the Lorentz structure analysis, the inverse of the quark propagator G(p) ^ and the scalar vertex T 
has such general forms. 


G(p) 1 = iy • pA(p^) + B(p^), 

(35) 

r™(p,0) = i-j ■pC(p^) + D(p^). 

(36) 

And according to the works of Ball and Chiu [50, 51], we take the following forms of Tj, 

,(q,p) and A^{q,p) 

rp(g,p) = Ea7p + (9+p)p[^7-(9 + p)Aa-*Ab] 

+(q + p)tj.[^i ■ (q + p)Ac - iAd], 


where 


2 

(37) 


and T = A, B, C, D. Here it should be pointed out that, there are also many other forms of vertex ansatz 
that can also satisfy the vector WTI. The main reasons that the authors of Ref. [46] chose the BC vertex is 
that: firstly, there is no kinematic singularity; secondly, the momenta p and q are symmetrical. The main 
purpose of their work is to discuss the differences of the results between the bare vertex and a corrected one. 

^Actually, the QCD vertex does not satisfy this identity, but rather the more complicated Slavnov-Taylor identity (STI), 
which can be derived from the Becchi-Rouet-Stora invariance of QCD and corresponds to the WTI of QED. The difficult 
challenge is to find an Ansatz for the unknown renormalised propagators, vertices, etc., of QCD which satisfy the DSEs and 
respect the STIs of the theory. Progress can be made by satisfying a subset of the DSEs and STIs, and then supplementing 
these with information gleaned from Lattice QCD, phenomenology, etc. However, without ghosts Eq. (33) is identical to the 
corresponding WTI of QED [83]. 
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Now substituting T^{q,p) into the DSE of the quark propagator Eq. (19), we can find the numerical 
solutions of the vector function A and the scalar function B of the quark propagator. Then by substituting 
T^{q,p), A^{q,p), A, and B into Eq. (32) we can solve the scalar vertex. This method can be generalized to 
the calculations of vacuum susceptibilities at finite temperature and finite quark chemical potential. In the 
following part we will show the details of the gluon propagators and the numerical results. 

3.2. Various vacuum susceptibilities 

In this part we discuss the calculations and results of various vacuum susceptibilities in detail, namely, 
the tensor, the vector, the axial-vector, the scalar, and the pseudo-scalar vacuum susceptibilities. 


3.2.1. The tensor vacuum susceptibility 

Studies of tensor vacuum susceptibility are closely related to the tensor charge of nucleon [10]. There 
are some calculations related to the tensor vacuum susceptibility, however, the results have shown that the 
theoretical treatments of this quantity is subtle and different approaches may lead to different outputs, even 
have different sign. Eor example, within the framework of QCD sum rules, H. X. He and X. D.Ji (1995, 1996) 
got about 0.002 GeV^ [11, 12], V. M. Belyaev and A. Oganesian (1997) got -0.008 GeV^ [13], L. S. Kisslinger 
(1999) got 0.0072 - 0.0127 GeV^ [15], A. P. Bakuleva and S. V. Mikhailov (2000) got -0.0055±0.0008 GeV^ 
for the non-local condensates sum rule and —0.0053 ± 0.0021 GeV^ for standard sum rule [16]; while W. 
Broniowski et al (1998) got -(0.0083 ~ 0.0104) GeV^ using the chiral constituent model [14], and H. T. Yang 
et al (2003) got -(0.0014 ^ 0.0020) GeV^ within global color symmetry model [17]. Clearly, in order to get a 
reliable theoretical prediction of the tensor charge, we should determine the tensor vacuum susceptibility as 
precise as possible. The work of Ref. [45] has provided a general treatment of the issue, and in Ref. [49] the 
calculations have been performed by employing the Ball-Chiu (BC) vertex [50, 51]. The results of Ref. [49] 
show that the tensor vacuum susceptibility obtained in the BC vertex approximation is reduced by about 10% 
compared to its rainbow-ladder approximation value, which means that the dressing effect of the quark-gluon 
vertex is not very large in the calculation of the tensor vacuum susceptibility within the DSEs framework. 

Through the previous general discussions of the vacuum susceptibilities (Sec. 2 of this paper) we know 
that, if we introduce a term that the quark tensor current coupled to a constant tensor external field 
q{x)(T^^q(x)Z^v{Z^i, is the tensor external field), using the QCD sum rules approach we can obtain the 
general expression of the tensor vacuum susceptibility (here is the chiral limit case, but it is easy to generalize 
to the case with nonzero current quark mass) 


^ {Tr[cr^<;Gr • ZG] - Tr[ar,cGP'^TP^ ■ ZG^^]} 
^ " Z,c{0\ : 9(0)<7(0) : JO) 

The expression in the earlier literatures is [14, 15, 16, 17] 

{Tr[ar,cGa ■ ZG] - Trla^cG^'^a ■ ZG^'^]} 
^ ~ Z,c(6l : g-(0)g(0) : JO) ’ 


(38) 


(39) 


Comparing these two equations it is easy to see that, the difference between and x'^ comes from the 
fact that in the previous calculations instead of using the exact dressed vertex and the perturbative one 
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r^J, they used the bare vertex in the course of which the non-perturbative and perturbative dressing 
effects on the vertex were ignored. Here we want to stress that, both the full and the perturbative tensor 
vertexes have complicated momentum dependence. Especially for the full tensor vertex, there may exist 
mass singularities of various meson bound states, so compared with the bare vertex it contains quite a lot 
of sophisticated non-perturbative effects of QCD. 

The discussions in Ref. [45] are under the rainbow-ladder approximation, and the renormalization pro¬ 
cedure is necessary. The renormalized version of the tensor vacuum susceptibility is (for the details please 
see Ref. [45]) 


^ _ {Tr[a^^GRTR ■ ZGr] - 
Xr — 


PT 

R 


■ ZG^r^]] 


(40) 

ZtZ„c( 0 | : 9fl(0)gfl(0) : 10 ) 

where Gr^ T/j, G^ and represent the renormalized versions of the full quark propagator, the full tensor 
vertex, the perturbative quark propagator and the perturbative tensor vertex, respectively. (Oj : qR(0)qR{Q) : 
jo) is the renormalized two-quark condensate. 

To show the numerical difference between s-nd we need to specify a model gluon propagator. The 
one Ref. [45] selected is the popular as well as famous Maris-Tandy model [87], which reads 


47r' 


D{k^) =—Dk'^e-^ -bdTT- 


7m TT 




2 In 


-f (l -b 




(41) 


here = [l — exp(—fc^/[4m(])] //c^, r = — 1, 7 ^ = 12/(33 — 2iVj) is the anomalous dimension of mass. 

For Nf = 4, = 0.234 GeV, and the renormalization point is chosen as /i = 19 GeV. As Ref. [87] 

pointed out, the parameters lu and rrit are not independent, to chose mt = 0.5 GeV and w = 0.3 GeV is to 
ensure that in the region 2 GeV^ < there is g^D{k^) ~ 47rQ;s(fc^), where as(fc^) is the “running” coupling 
constant of QCD. In Ref. [87] they chose D = 1.25 GeV^, and there studies show that this set of parameters 
can give very good results for the properties of pion and kaon. 

Once the model gluon propagator is chosen, we can then calculate the quark propagator and the vertex 
through numerical iteration, as shown in Fig. 1 and Fig. 2 respectively. The final numerical results of the 
tensor vacuum susceptibilities are 


X^a = -0.0027 GeV^ , y'^a = -0.0031 GeV^ 


(42) 


here a = -4(6] : qR{0)qR{0) : | 0 ). 

We see that, and x'^ are indeed different, but the difference is not very large. Let us now analyze the 
reasons for this. As stated above, the difference between them comes from the choices of vertexes. We know 
that in the perturbative region the bare vertex can be regarded as the zeroth approximation of the dressed 
vertex T^j/, and therefore it is physically interesting to analyze the effects caused by their differences in the 
non-perturbative region. It can be seen clearly from Fig. 2 that, for large enough momentum square s = 

Ai(s,/i^) tends to 1 while A 2 (s,/r^) tends to 0 , which means that in the large momentum region tends 
to however, for small s, Ai(s, p^) differs from 1 and differs from 0, although the differences 

are not significant. Therefore, the dressed vertex rpi^(p, 0) obtained in the rainbow-ladder approximation of 
the DSEs approach is similar to that of the bare vertex in the small momentum region, and it is just this 
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Figure 1: The scalar functions yl(s, /r^) and B{s, [GeV] of the dressed quark propagator, taken from Ref. [45]. Here s = p^. 




Figure 2: The tensor vertex functions Ai(s,p?) and A 2 (s,/r^)[GeV ^], taken from Ref. [45]. Here s = ■ 

small difference that accounts for the small numerical difference between the numerical values of s-nd x'^ ■ 
The small magnitude of the numerical difference between x^ x'^ suggests that the non-perturbative 
dressing effects on the dressed tensor vertex are not important, and therefore in this case to use as 
a replacement of 0) is a rather good approximation. This is quite different from the case of vector 

vacuum susceptibility [40]. 

3.2.2. The vector and axial-vector vacuum susceptibilities 

For the vector vacuum susceptibility xv and the axial-vector vacuum susceptibility xa, in order to 
distinguish them from the later vertex wave function Xu, here we use different symbols Ky and ka instead. 
According to the previous derivations, ny can also be written as 

Ky{qq) = Jd^zJ^ [n^^iP) - n^’^iP)] , (43) 
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where is the vector vacuum polarization 

pA 

ul,{p) = -Z2 tr[j^xAq;P)], (44) 

Jq 

while for 

=jdAj^ ^,p. Mz)Vs^} _ „A.PT^p^] _ ,45) 

with the axial-vector vacuum polarization 

pA 

^%{P) = -Z2l tr[x^XtiXbu{q]P)], (46) 

Jq 

in which the trace is to be taken in color and Dirac space. V^{z) and V 5 ^{z) represent the variable vector and 
axial-vector external fields, the exact vector and axial-vector vertex wave functions, which can 

be written as C?rG using the no external field versions of the full quark propagator G and the corresponding 
vertex functions F. Z 2 is the renormalization constant of the quark wave function. According to the WTI, 
the vertex renormalization constant is also equal to Z 2 . {qq) is he chiral quark condensate. 

As can be seen from Eqs. (43) and (45), in order to obtain the non-perturbative vector and axial-vector 
vacuum susceptibilities, we should subtract li^P'^{P) and FI^^^^(P) respectively, which arises from the 
perturbative effects. In other words, when calculating the vacuum susceptibilities, the mean values of the 
perturbative vacuum should be subtracted. This is very important and is deserving of additional attention. 
It is well known that the separation of the perturbative and the non-perturbative contributions from the 
mean values of the vacuum is somewhat arbitrary. Usually, this arbitrariness is avoided by introducing some 
normalization point [123]. In such a formula, the condensates will depend on the choice of the normalization 
point. There are also other methods besides this one, for example, in studies of the mixed quark-gluon 
condensate, the authors of Ref. [124] identified the perturbative vacuum with the Wigner vacuum, since both 
of them are trivial in the sense that there are no chiral symmetry breaking and confinement, in contrast to 
Nambu-Goldstone vacuum (the non-trivial vacuum) which corresponds to DCSB (more details can be found 
in Ref. [124]). In Refs. [40, 44, 45], the authors adopted the viewpoint of Ref. [124] to calculate the vector, 
the axial-vector, and the tensor vacuum susceptibilities in the framework of rainbow-ladder approximation of 
the DSEs approach. For example, in the calculation of the vector vacuum susceptibility, in Eq. (43) 

is calculated in the Nambu-Goldstone vacuum configuration, while 11^^^^(P) is calculated in the Wigner 
vacuum configuration. It is obvious that this calculation depends on the rainbow-ladder approximation of 
the DSEs approach. In the literature, there are few theoretical studies related to the vector and axial-vector 
vacuum susceptibilities, among them, L. S. Kisslinger determined them using a three-point formalism within 
the method of QGD sum rules [15], M. Harada et al discussed the effective degrees of freedom at chiral 
restoration and the vector manifestation in hidden local symmetry theory [125], and K. Jo et al calculated 
vector susceptibility and QCD phase transition in anti-de Sitter (AdS)/QGD models [126]. In the following, 
we will show how the authors of Ref. [48] obtain model independent results for the vector and axial-vector 
vacuum susceptibilities using the vector and axial-vector WTIs. 

Let us discuss the vector vacuum susceptibility Ky first. The conservation of the vector current ensures 
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that the vector vacuum polarization is purely transverse, namely, 




p2 


If we take the constant external field limit V^(z) = V^, and use the following integration equation 


dH S\l) Pd^f^k,pJ) = \j dH 6\l)p-kf{k,p,l), 


we can get 


Kv{qq) = I [n^(P' = 0) - = 0) 


( 47 ) 


(48) 


(49) 


It can be seen that the vector vacuum susceptibility is closely related to the vector vacuum polarization at 
zero total momentum. Now contracting both sides of Eq. (44) with we then have 

rA 


n^(p2=0) = -^y tr [ 7 ^x^(q;P = 0 )] . 


(50) 


The case for 11^= 0) is similar. 

In the chiral limit case, the wave function of the vector vertex satisfy the WTI 


P) = G{q-) - G{q+), 


(51) 


where q± = q± P/2. Expanding the right side of Eq. (51) to Pf^ and taking the limit P^ —>■ 0 leads to the 
Ward identity 


Xf,{q\P = 0) = i 


dG{q) 

dqu 


(52) 


Substituting Eq. (52) into the right side of Eq. (50), and adopting the following parametrization of the quark 
propagator G{q) 


G{q) = 


1 


ix ■ kA{q^) + B{q^) 


= -ij- kav{q^) + CTs(q^), 


(53) 


we obtain 


n^(p2 = 0) = tr 


. dG{q) 




dqt^ 

= -8Z,J [2av{q^) + q^^^0^] 


(54) 


,q dq^ 

Since for large q^, cry ~ q~^, this integral is quadratically divergent. However, this divergence is not 
genuine. Note that the integrand is a total divergence, so the above integral vanishes if a translation 
invariant regularization is adopted,® and hence 


n^(p^ = 0) = 0. 


(55) 


Eor similar reasons we can get H^’'^^(P^ = 0) = 0 too. Thus we may draw the conclusion that the vector 
vacuum susceptibility is zero so long as the Ward identity is satisfied. In other words, the vanishing of the 


® Actually, in numerical calculations of the vacuum polarization using DSEs people usually employ a cutoff to regularize the 
ultraviolet divergence in Eq. (44), for the details please see Ref. [48] and references therein. 
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vector vacuum susceptibility is a direct consequence of gauge invariance. The authors of Ref. [40] also made 
use of the vector Ward identity in constructing the vertex function, but at that time they did not realize that 
the vector vacuum susceptibility can be expressed with such a simple formula, but their numerical result (of 
the order 10“^) is very close to this model independent one. 

The case for the axial-vector vacuum polarization is a little bit more complicated, since the axial-vector 
vertex contains a massless bound state pole. Nevertheless, the analysis is still straightforward, and the 
general form of the axial-vector vacuum polarization can be written as 

( 56 ) 

which contains both a transverse part and a longitudinal part. Following the same procedure of deducing 
the vector vacuum susceptibility, we can get 


KA{qq) = J [n^(p2 = 0) - = 0) 


nf (p2 = 0) - = 0) 


(57) 


From Eqs. (56) and (46), we know that 

n^(P') - ^nf(p2) = -Z2 j\r[j5l^X5AQ-,P)] • (58) 

Now we will determine the axial-vector vacuum polarization at P^ = 0. In the chiral limit, the WTI for 
the axial-vector vertex can be expressed as 


-iP^iXbM^P) = G('7-e)75 + 75^(9-) • (59) 

where we also express the identity in terms of the dressed axial-vector vertex wave function Xbfi instead of 
the dressed vertex Fs^. Since Xb^^ possesses a longitudinal massless bound state pole, its generally expression 
is then [54] 

X5m(9; p) = 75x^(9; p) + p) + p) ■ (60) 

where x^(g; P) = ^fj.FR+x-qqfj,GR-a^^q^HR, functions Fr, Gr, Hr are regular at P^ 0, P/xXs/i ~ O(P^), 
/,r is the pion decay constant in the chiral limit, and XTr{q\P) is the canonically normalized Bethe-Salpeter 
wave functions of the massless bound state which take a general form 


Xniq',P) = 275 [iE^{q;P) -ky- PP^(g;P) + x-qq-PG^{q; P) + a^i,qf,P^H^{q-, P)]. 
In the chiral limit, /tt is defined as 

f^P^ = Z 2 [ tr [757A«X7r(g; T")] • 

Jq 

It is apparent that as P —>■ 0, we can expand the right side of Eq. (62) to O(P^), and then 

/^ = -24Z2^ (^P,(g;0) + ir72G^(g;0)^ . 

Contracting both sides of Eq. (58) with the projector — 4 " , we can get 


3n^(p2) -t snf (p2) = -Z 2 / tr 

J a 


PUibiriinbXv^F P) + Xhv{q\ P) + ^^^xAq\ P)) 


(61) 


(62) 


(63) 


(64) 
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Not let us focus on the limit of Eq. (64) when —>■ 0. After some algebra we find that, the first term in the 
right side of Eq. (64) vanishes, and the second term also vanishes thanks to ~ 0{P^), while the third 

term is found to be —24/^A^c-^2 + ^q^GT^(q; 0)), which equals 3/J after making use of Eq. (63). 

Therefore we obtain 

n^(p2 = o) + nf(p2 = o) = /2. (65) 


Similarly, contracting both sides of Eq. (58) with the projector — C (C 7 ^ 4) gives 


m^{p^) + {C-mt{P^) = -Z 2 ttr 

J Q 


,Rt„. u\ I u\ I 

p2 


Pp.ui5i^.{ibXu (<?; P) + Xbu{q] P) 


x^{q\P)) 


( 66 ) 


In the limit —>• 0, the first term of Eq. (66) is (4 — C)ANcZ 2 {FR{q', 0) + jq^Gn{q; 0)), the second term 

does not contribute, while the third term is {( — 1)/^. Now using Eq. (65) we then get 


n^(p2 = 0) = AN,Z2 j [Efl(g;0) + \q^GR{qM- 


(67) 


Note that the result is independent of (^. 

Now substituting the general form of Xb^i in Eq. (60) into Eq. (59), and expand both sides to 0(1) and 
0{Pfj) respectively, we can get the corresponding Goldberger-Treiman relation. Eunctions Pij(g; 0), Gfi^q-, 0) 
can be expressed using the vector part of the quark propagator and the BS wave function of pion (details 
can be found in Ref. [54]) 


Pfl(g;0) = -ay(g2)-2/,P^(g;0), 

Gfl(g;0) = -2^^^^-2/^G^(g;0). (68) 

Substituting Eq. (68) into Eq. (67), and using Eq. (63), we get 

n^(p2 = 0) = -2N,Z2 ^ [2av{q^) + + /' = fl (69) 

To obtain this result we have used the previous discussion that when using a translation invariant regular¬ 
ization, the integral to q in Eq. (69) is 0. 

Combining Eqs. (65) and (69), we have 

n^(p2 = o) = /2, n^(p2 = o) = o. (70) 

The subtraction term (P^ = 0) and n^’^^(P^ = 0) can be obtained using similar methods. Note that 

the perturbative axial-vector vertex function has no pion pole term. We obtain 

nAPT(p2 ^ 0) = 0, n)^’-^^(p2 = 0) = 0. (7i) 


The axial-vector vacuum susceptibility is then 

i^A{qq) = \fl- 

Erom Eqs. (55) and (70), we found that in the chiral limit case 

n^(p2 = o)-n^((p2 = o) = -/2. 


(72) 


(73) 
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which is just the result from Weinberg sum rules [127]. 

In summary, in this part we have derived model independent results of the vector and the axial-vector 
vacuum susceptibilities, which will play an important role in the related calculations of the QCD sum rules 
external field method. Using the vector and the axial-vector WTIs, Ref. [48] has demonstrated that in the 
chiral limit the vector vacuum susceptibility is 0, while the axial-vector vacuum susceptibility equals 3/4 of 
the square of the pion decay constant, which is the same as the result from Weinberg sum rules. 


3.2.3. The scalar vacuum susceptibility 

Studies on the QCD phase transitions play important roles in the studies of the evolution of the early 
universe and the high energy heavy ion collisions. The results of the Lattice QCD show that the QCD “phase 
transition” in the early universe is not a real phase transition, but a smooth crossover [58] . In Ref. [58] they 
adopted the scalar vacuum susceptibility x an indicator, and studied the temperature dependence of y. 
The scalar vacuum susceptibility is very important in studies of QCD phase transitions, however, it is known 
that there is a square divergence in the x(T) calculation. In order to eliminate it, a usual treatment is to 
subtract x(T = 0) from x(T), since the divergence only exists in y(T = 0), and the temperature dependent 
part is not divergent [57, 58]. Obviously, such method is not applicable in the zero temperature case, but 
the scalar vacuum susceptibility is meaningful at zero temperature. We will give a new treatment to the 
square divergence in x(T) later. 

In QCD, the two-quark condensate associated with / flavor of quarks is defined as [128] 

rA 


{qq)f{mf] /i. A) = AjWtro / G/(g; mf, p) , 

Jn 


(74) 


where mf{fi) is the renormalized current quark mass, with /i the renormalization scale; Z 4 is the renor¬ 
malization constant of the mass term in the Lagrangian, which depends implicitly on the gauge parameter. 
In the chiral limit case the two-quark condensate can be used as an order parameter to describe the chiral 
phase transition, the scalar vacuum susceptibility (or chiral susceptibility) just characterize its response to 
the change of quark mass [31, 129] 


x(m) := 


d 


dm[p) 


{qq){m;p.,A) 


7h—0 


d 

= Za[pl, A)A^ctrD J "i; m) 


(75) 


rh=0 


here m is the auxiliary current quark mass that can be regarded as renormalization invariant. The order of 
integration and differentiation can be interchanged because the theory is properly regularized. Now use the 
Ward identity 


then 


—G(g; m; p.) = -G{q; m; n)T{q, 0; /r)G(q; m; p) 


XM = -Ziin, A)iVctrD f G{q] 0; /i)r(g, 0; /r)G(g; 0; fi), 

J n 


where T is the renormalized fully-dressed scalar vertex, which satisfies the inhomogeneous BSE 

rA 


r(fc, P;/r) = ZaId 


[S{q+)T^{q, P)S{q.)UKi:{q, fc; P). 


(77) 


(78) 
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with k the relative and P the total momenta of the quark-anti-quark pair; q± = q ± P/2; r, s, t, u represent 
color and Dirac indices; and K is the fully-amputated quark-anti-quark scattering matrix. 

The quantity m^no(P) is renormalization point invariant in QCD, where no(P) is the scalar vacuum 
polarization 

noiP;fi) = ZiN,tTD [ G(g+)r(g,P)G(g-). (79) 

J q 

To compare with Eq. (77), we can get the general result 

X(p) =-no(P = 0 ;to = 0 ,^). (80) 

Hitherto we have not specified a regularization procedure for the scalar vacuum susceptibility. In this 
connection it is noteworthy that if a hard cutoff is used, then in the chiral limit of a non-interacting theory 

-nO(P = 0;m = 0) = ^A2. (81) 

Ref. [128] pointed out that, this result can be traced to the dependence on current quark mass in Eq. (74). 
On the other hand, the Pauli-Villars regularization would yield zero, which is then the procedure that the 
authors of Ref. [46] recommend and employ in models that preserve the one-loop renormalization-group 
behavior of QCD, and in this case the scalar vacuum susceptibility is defined self-consistent. 

For the calculation of xik) only need the scalar vertex at P = 0, at which total momentum it has the 
general form 

ro(fc,0; m) = ij ■ kC{k^]m) + D{k‘^]m). (82) 

Owing to the Ward identity, Eq. (76), it is not necessary to solve the inhomogeneous BSE since 

G(fc^;m) = ——H(fc^;m); D{k^\m) = — — B(k^-,m). (83) 

dm am 

Hence a solution of the gap equation suffices completely to fix ro(fc,0). 

A simplified form of the renormalization group improved effective interaction is employed in Ref. [46]; 

viz., they retain only that piece which expresses the long-range behavior [85, 87] 

2 2 

g'^Df.^iq) = ATr'^t^^{q)D^ exp(-^), (84) 

with tf^v{q) the transverse momentum projection operator. Eq. (84) delivers an ultraviolet finite model gap 
equation. Hence, the regularization mass scale A can be moved to infinity and the renormalization constants 
equal to one. 

The active parameters in Eq. (84) are D and to, but they are not independent. In reconsidering a 
renormalization group improved rainbow-ladder fit to a selection of ground state observables [87], Ref. [93] 
noted that a change in D is compensated by an alteration of w. This feature has further been elucidated and 
exploited in Refs. [130, 131]. For the range [0.3,0.5] GeV, the fitted low-energy observables are approximately 
constant along the trajectory 

wP> = (0.8GeV)3 =: (85) 

Herein, we employ w = 0.5 GeV, and therefore D = m^fuj = 1.0 GeV^ corresponds to what might be called 
the real-world reference value for the bare vertex. Here we should note that, if we choose Eq. (84) as the 
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model gluon propagator, but for the quark-gluon vertex using the BC ansatz Eq. (37), there is no reason 
to suppose that w and D still satisfy the relation in Eq. (85), since the BC vertex contains considerable 
unknown non-perturbative information of interaction. Therefore, we need to refit the parameter D. After 
some tries, it has been found that if one assumes D = 0.5 GeV^, we can get the value of the two-quark 
condensate shown in the fourth column of Table 1. Hence we postulate that D = 0.5 GeV^ is the real world 
reference value under the BC vertex. 

Table 1: Results obtained for selected quantities with u) = 0.5 GeV, and the vertex and D parameter value indicated: A(0), 
M(0) are p = 0 in-vacuum values of the scalar functions A{p^) and M(p^) = the quark condensate is defined 

with m = 0 in Eq. (74); and x i® obtained from Eq. (90). A(0) is dimensionless, but all other entries are quoted in GeV. The 
calculations reported herein were performed in the chiral limit, and all are taken from Ref. [46] 


Vertex 

Vd 

A(0) 

M(0) 

-((q-q)0)V3 

Vx 

Rainbow-Ladder 

1 

1.3 

0.40 

0.25 

0.39 

Ball-Chin 

1/V2 

1.1 

0.28 

0.26 

0.28 


As discussed above, Eq. (81) demonstrates that a regularization procedure is needed for the vacuum 
polarization, even if dealing with a free field theory. For QCD, in which Schwinger functions carry anomalous 
dimensions, the Pauli-Villars scheme works well. However, this procedure fails when the vacuum polarization 
is evaluated with Schwinger functions generated by Eq. (84), where the interaction does not preserve the 
one-loop behavior of QCD, but delivers an ultraviolet finite gap equation instead. In the following we will 
construct a simple alternative. 

The vector vacuum polarization is given by 

pA 

= ^cHd / i7A.G(g+)ir,(g,P)G(g_), (86) 

J q 

here T^ is the vector quark-anti-quark vertex that satisfies an inhomogeneous BSE similar to Eq. (78). This 
polarization couples to the photon and hence must be transverse while contain no mass term. Consider 
therefore [48] 

1 N 

■^^lf,iP = 0) = ^tro J ijf,G{q)iT^{q,0)G{q) 

= (87) 

where the vector Ward identity was used in the second line and we have written the dressed quark propagator 
in the following form 

G{p) = -i"f ■ pavip^) + crsip^). 

If we consider a free field and use a hard cutoff, then Eq. (87) becomes 


( 88 ) 






namely, the same sort of divergence is encountered as in the scalar vacuum polarization, Eq. (81). Naturally, 
any regularization scheme that preserves the vector Ward-Takahashi would yield zero as the result [132]. So, 
considering Eq. (87), we can define the regularized scalar vacuum susceptibility as 

X = -no(0; m = 0, A) + inJC^(0; m = 0, A) (90) 

N 1'°° ( 

= J dss |a’(s)[o-s(s)^ - scrs(s)^] 

+2s C'(s)crv(s)o’s(s) + 2cry(s) + scry(s)|. 

Here, as long as each term in Eq. (90) is regularized independently in a valid fashion, then this is a “nugatory” 
transformation; since the photon is massless and hence in any valid scheme there must be n^^(0; 0, A) = 0. 
In other words, with this definition we can move on to calculate the susceptibility without specifying a 
particular regularization scheme. The last two contributions in the integrand act as the Pauli-Villars terms, 
and Z 4 = 1 in the ultraviolet-finite model. 

Now we will discuss the numerical results. Owing to the Gaussian form of the gluon propagator, Eq. (84), 
all the relevant integrations converge rapidly. Fig. 3 shows the functions obtained by solving the gap equation, 
and Fig. 4 gives the results which describe the P = 0 scalar vertex. All the results were obtained with the 
appropriate real world values of the interaction strength: D = 1 GeV^ for the bare vertex, and D = 0.5 GeV^ 
for the BG vertex, as listed in Table 1. 



Figure 3: The dressed quark propagator, taken from Ref. [46]. Left panel: A{p^) = 1/Z(jp), where Z{p^) is the wave function 
renormalization function; Right panel: B{p^) [GeV], the scalar piece of the dressed-quark self-energy. In both panels, the 
dashed curve is calculated in the rainbow-ladder truncation, with X = DIhP = 4; while the solid curve is calculated with the 
BC vertex ansatz, and X = 2. 

It can be seen clearly from Fig. 3 that, the BC vertex ansatz has a quantitative impact on the magnitude 
and pointwise evolution of the gap equations solution. This can be anticipated from Ref. [133]. Moreover, 
the pattern of behavior can be understood from Ref. [112]: the feedback arising through the term in 
the BC vertex always acts to alter the domain upon which AijP') and Bijp') differ significantly in magnitude 
from their respective free-particle values, especially in the intermediate momentum region. The behaviors of 
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Figure 4: P = 0 scalar vertex, taken from Ref. [46]. Left panel: C'(p^); Right panel: D{p^). In both panels, the dashed curve is 
calculated in the rainbow-ladder truncation, with X = D juP" = 4; while the solid curve is calculated with the BC vertex ansatz, 
and X = 2. 

C{p^) and D{p^) shown in Fig. 4 can be understood easily from those of A{p^) and B{p^), so that there is 
no need to discuss it further. 

The integrand in Eq. (90) is depicted in Fig. 5 for each vertex ansatz at the associated real world 
interaction strength. The resulting scalar vacuum susceptibilities are presented in the fifth column of Table 1. 
It can be found from this figure that the regularization has served to eliminate the far-ultraviolet tail of the 
integrand, thereby ensuring convergence of the integral. We have varied the detailed form of the regularizing 
subtraction, namely, using free-held propagators and vertices instead of gap and BSE solutions, and the 
behavior of the integrand changes a little. Moreover, at real world values of the interaction strength for 
both ansatzs the integrands have negative support in the infrared and positive support for > (O.SGeE)^. 
These results tell us that when studying the scalar vacuum susceptibility, to go beyond the rainbow-ladder 
truncation is important and necessary, which then deserves further investigation. 

In Fig. 6 we depict the evolution of the scalar vacuum susceptibility with increasing interaction strength, 
I = D/uf^. For 1 = 0 there is no interaction and then the “vacuum” is unperturbed by a small change in 
the current quark mass, hence the susceptibility keeps zero. With increasing X the susceptibility will grow 
since the interaction in the qq channel is attractive and therefore magnifies the associated pairing. This is 
equivalent to stating that the scalar vertex is enhanced above its free field value. Then the growth continues 
and accelerates until at some critical value X^ the susceptibility becomes infinite, namely, a divergence appear. 
The critical values are X^. = 1.93 for the rainbow-ladder approximation, and X,, = 1.41 for the BC vertex 
ansatz. 

This divergence can be understood in this way: the models we have defined contain a dimensionless 
parameter I, which characterizes the interaction strength, and a current quark mass, which is an explicit 
source of chiral symmetry breaking. In the general theory of phase transitions, the latter is analogous to an 
external magnetic field while 1/X is kindred to a temperature. Consider the free energy for such theories 
f{t, m), where t = [1—Xc/X]. If such a theory possesses a second-order phase transition, then the free energy 
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(GeV^) 


Figure 5: Integrand in Eq. (90), taken from Ref. [46]. Dashed curve: calculated in rainbow-ladder truncation, with X = DjuP' = 
4; solid curve: calculated with BC vertex ansatz, and X = 2. 

is a homogeneous function of its arguments in the neighborhood of t = 0 = m. From this it follows that the 
theory’s magnetization exhibits the following behavior (e.g., see the Appendix of Ref. [27] for the details): 

M{t,Q) = t^, (91) 

and the associated magnetic susceptibility evolves in this way 

M(0,m) oc (92) 

For a mean-field theory, /3 = 1/2, and (5 = 3. 

Now the nature of the critical interaction strength is easy to understand. In the class of theories we 
are considering, the quark condensate is analogous to the magnetization, and it is attended by the scalar 
vacuum susceptibility. For X < Ic, the interaction is insufficient to generate a nonzero scalar term in the 
dressed-quark self-energy in the absence of a current quark mass; namely, DCSB is impossible and the model 
realizes chiral symmetry in the Wigner mode. But the situation changes at Ic, and for I > Ic a, B ^ 0 
solution is always possible. Moreover, the behaviors of the susceptibility show that each model undergoes 
a second-order phase transition and realizes chiral symmetry in the Nambu-Goldstone mode for interaction 
strengths above their respective values of Ic- These observations emphasize the usefulness of the scalar 
vacuum susceptibility. With the bare vertex people can construct the pressure explicitly,® and thereby show 

®The pressure is defineci as the negative of the effective action, and then the ground state of a system is that configuration for 
which the pressure is a global maximum. For the rainbow-ladder approximation, the corresponding pressure is just constructed 
from the famous CJT effective action [134] . 
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Figure 6: Dependence of the scalar vacuum susceptibility on the interaction strength in Eq. (84), taken from Ref. [46]; viz., 
X ■.= DjuP". dashed curve is the RL vertex, while solid curve is BC vertex 

that for X>Xc the DCSB solution is dynamically favored because it corresponds to the configuration of the 
maximum pressure. Nevertheless, the diagrammatic content of the BC vertex ansatz is not knowable and 
hence an expression for the pressure cannot be derived. In this case we may rely on the behavior of some 
susceptibilities to conclude that DCSB is favored, as illustrated within many model studies. Once past the 
critical point, the susceptibility decreases as X increases. This is because the magnitude of the condensate 
order parameter grows in tandem with X and, therefore, the influence of any perturbation associated with a 
current quark mass must steadily diminish. We will do some related discussions in Sec. 4 of this paper. For 
the theoretical studies of the scalar and the following pseudo-scalar vacuum susceptibilities using the QCD 
sum rules, please refer to Ref. [15]. 

3 . 2 . 4 . The pseudo-scalar vacuum susceptibility 

The color-singlet current-current correlators (or equivalently, the associated vacuum polarizations) are 
directly related to physical observables, and hence play important roles in QCD. For example, the vector vac¬ 
uum polarization is coupled to real and virtual photons, so that it is basic to the analysis and understanding 
of the process e+e“ —>■ hadrons [135, 136]. The standard Lattice QCD methodology related to the correlators 
can be found in Ref. [137]; on the other hand, correlators are also amenable to analysis via the operator 
product expansion and are therefore fundamental in the application of QCD sum rules. In the latter connec¬ 
tion, the pseudo-scalar vacuum susceptibility (or the pion susceptibility) plays a role in the QCD sum rules 
estimate of numerous meson-hadron couplings, for example, the strong and parity-violating pion-nucleon 
couplings, Pttnn and respectively [7, 138, 139]. Furthermore, the pseudo-scalar vacuum susceptibility 

is also a probe of QCD vacuum structure as the scalar susceptibility which is discussed in Sec. 3.2.3 of this 
paper [46], while its analysis is more subtle, with conflicts and misconceptions being common in different 
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model calculations [7, 8 , 33, 139, 140]. 

In this part we derive a model independent result for the pseudo-scalar vacuum susceptibility using the 
isovector-pseudo-scalar vacuum polarization [47] , which can be written as 

nf (P;C) = Z4iV.tr ^^^75r^G(<7+)*rg(Q;P)G(g_). (93) 

where ^ represents the renormalization scale, the trace is over flavor and spinor indices, is the Pauli matri¬ 
ces of the isospin, P is the total momentum of the quark-anti-quark pair; G{q±) represents the dressed-quark 
propagator, and Ps is the fully dressed pseudo-scalar vertex, both of which depend on the renormalization 
point. Physical quantities obtained from Eq. (93) are then gauge invariant. The propagator can be obtained 
from the gap equation, Eq. (19). 

The pseudo-scalar vertex is solved from an inhomogeneous BSE, namely 

[ri(fc;P)]t„ = l\xi(.q-,P)]srKi:{q,t,P), (94) 

here k is the relative momentum, r,s,t^ and u represent color, flavor, and spinor indices, and 


X^(fc;P)=G(fc+)r^(fc;P)G(fc_), 


(95) 


with k± = fc±P/2, and K{q, k] P) is the fully amputated two-particle irreducible quark-anti-quark scattering 
kernel. Now, let us consider the case with a space-time independent pseudo-scalar source, S 5 4 associated 
with the term 

J d‘^xq{x)^-/5TS5q{x) ( 96 ) 

in the action, we can define a pseudo-scalar vacuum condensate, the gauge-invariant as well as properly 
renormalized form of which in QCD is 

i 

{qq)^{s5,m-,(,A) = Z4Ncti J -75r'"G(g; S5, m; C) • ( 97 ) 

We can see that it is analogous to the vacuum quark condensate, Eq. (74). When m = 0, the DCSB can be 
expressed via the Higgs mechanism as 


- {qq)c = lim (gg)(m;C, A) 4 0. 

^ m—>-0 


(98) 


Eq. (98) then defines what we mean by an isoscalar-scalar configuration: isovector-pseudo-scalar correlations 
are by convention measured with respect to this configuration. These observations highlight the importance 
of the pseudo-scalar vacuum susceptibility 


xr(C) = ^( 99 ) 5 ( 55 ,w;C, A) 


S5=0 


Then following the previous discussions, we can know that 


(99) 


xr(C) = -2nr(H = 0;s5 = 0,m,C). 


( 100 ) 


Hitherto we have not specified a regularization procedure for the susceptibility, actually it can rigorously be 
defined via a Pauli-Villars procedure, as discussed in Ref. [46]. 
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Now let us discuss the value of in ths neighborhood of the chiral limit, therein we may write [54] 

0) = 0) + o). (loi) 

where E^(k-,P) is a part of the inhomogeneous pseudo-scalar vertex that is regular as + nri^ —>■ 0 , 
TT^fj,(k;P) is the pion bound state’s canonically normalized Bethe-Salpeter amplitude, and r-wiC) determined 

by ^ 

= {^\Q^l5T''q\Trf^) = ZiN^tr -P) (102) 

is the residue of this bound state in the inhomogeneous pseudo-scalar vertex. Using this, the pion’s leptonic 
decay constant can be expressed as 

S^''f^Pcr = {0\q^'l5laT''q\7r>^) = Z 2 ]pblaT''Xn{q'.P) ■ ( 103 ) 

E^{k]P) can be determine via the axial-vector WTI 

P.r^,(fc; P) + 2m(C)^^^(fc; P) = G{k+)-^\ix,T^^ + \iX5T^G{k .)-^, (104) 

with Fg^ the inhomogeneous axial-vector vertex. At P = 0 with m 7 ^ 0 there is no pole contribution on the 
left side and hence 

m(C)Pf(fc;P = 0) = P(fc^m;C^). (105) 

in other words, this regular piece of the pseudo-scalar vertex is completely determined by the scalar part 
of the m 7 ^ 0 quark self-energy. Using the systematic, non-perturbative, and symmetry-preserving DSE 
truncation scheme introduced in Refs. [84, 141], we can verify this equation order-by-order via the gap and 
Bethe-Salpeter equations. 

Now substituting Eq. (101) into Eq. (100), we then get 

Xf ^'^'^X5(C) (106) 

X 5 (C) = xl{C) + xf{0+O{m), (107) 


so that, in the neighbourhood of m = 0 , the pseudo-scalar susceptibility splits into a sum of two terms: the 
first one expresses the contribution of the pion pole, 0 (m“^), and can be expressed in the closed form 

2r^(C)^ m=o {qq)l 


xUO = 2 

of which the last equality is proved in Ref. [54]; the second term 0{rh^) can be determined via 


mi. m(C) 


m{0x5(CW' ”=° ZiNctr J G{q)^X5T'" B{q^ ,m)G{q) = (5'"‘'(gg)(m; C, A), 

where {qq){m\C,, h) is the vacuum quark condensate defined in Eq. (74), and this entails 

xiiCim) = x(C) -bO(m), 


(108) 


(109) 


( 110 ) 


where x(C) is the scalar vacuum susceptibility defined in Eq. (75). Hence we arrive at a model independent 
consequence of chiral symmetry and the pattern by which it is broken in QCD, namely. 


{qq)c 


X5(C) = - 


m(0 


+ xiO + Oim). 


(Ill) 
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Table 2: Pseudo-scalar vacuum susceptibility and related quantities computed using the two kernels of the BSE described in 
connection with Eqs. (37) and (84), also similar to those of Table 1. Dimensioned quantities are listed in GeV, and all are taken 
from Ref. [46] 


Vertex 

Vd 

UJ 

-((g-g)0)i/3 

fO 

J TT 

m 


vV 

Rainbow-Ladder 

1 

0.5 

0.25 

0.091 

0.0050 

1.77 

0.39 

Ball-Chiu 

1/V2 

0.5 

0.26 

0.11 

0.0064 

1.66 

0.28 


Now we show the numerical results of Ref. [47] in Table 2, which were computed from two models for the 
gap equations kernel. 

Eq. (Ill) is a meaningful result. Recall that, in the absence of a current-quark mass, the two-flavor 
action has a SUl(2) SU_r( 2) symmetry and, moreover, that ascribing scalar-isoscalar quantum numbers 
to the QCD vacuum is a convention contingent upon the form of the current-quark mass term. It follows 
that, the massless action cannot distinguish between the continuum of sources specified by 


constant x 



\ 8 \ e [0,2^). 


( 112 ) 


Therefore, the regular part of the vacuum susceptibility must be identical when measured as the response 
to any one of these sources, so that there is XR = X foi' choices of 9. This is the content of the so-called 
“Mexican hat” potential, which is used in building effective models for QCD. The magnitude of x depends 
on whether the chiral symmetry is dynamically broken and the strength of the interaction as measured with 
respect to the critical value required for DCSB, as discussed in Sec. 3.2.3 of this paper [46]. When the 
symmetry is dynamically broken, then the Goldstone modes appear, by convention, in the pseudo-scalar- 
isovector channel, and thus the pole contributions appear in xs but not in the chiral susceptibility. It is 
valid to draw an analogy with the Weinberg sum rules [48, 127]. Eq. (Ill) also provides a novel, model 
independent perspective on a mismatch between the evaluation of the pseudo-scalar vacuum susceptibility 
using either a two-point or a three-point sum rules. The two-point study in Ref. [140] produces the pion pole 
contribution, Xg, which is also the piece emphasized in Ref. [33], whereas a three-point method in Ref. [139] 
isolates the regular piece, x|^, since a vacuum saturation ansatz is implemented in the derivation. Thus their 
analyses are not essentially in conflict, but emphasize different and independent pieces of the susceptibility, 
which can be distinguished. However, only the regular piece should be retained in a sum rules estimate 
of the pion-nucleon coupling constants [138]. We note in closing that other vacuum susceptibilities can be 
analyzed similarly. 


4. Susceptibilities and the chiral phase transition of QCD 

The chiral phase transition of QCD and a possible CEP on the QCD phase diagram have been drawing 
much attention on both experimental and theoretical sides; for example, beam energy scan (BES) program 
at the RHIC [142, 143, 144, 145], and some theoretical studies [23, 57, 58, 60, 146, 147, 148, 149, 150, 
151, 152, 153, 154, 155, 156, 157, 158, 159]. To study the chiral phase transition of QCD, in the chiral 
limit the two-quark condensate can serve as an order parameter. Nevertheless, the current quark mass is 
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small but not zero in the real world, so when studying the partial restoration of chiral symmetry at finite 
temperature and quark chemical potential, people also need to resort to some alternative order parameters. 
When there is a coexistence of different phases, one of the candidates is the pressure difference (the bag 
constant) [160, 161, 162, 163]. As discussed above, various susceptibilities are the parameters that related 
to the linear responses of the system to the external helds, which can then characterize the properties of the 
system, as well as to serve as indicators of corresponding phase transitions. In this section, we focus on the 
applications of some QCD susceptibilities to the chiral phase transition of QCD. 

In quantum field theory, the dynamic properties of a system are fully characterized by the generating 
functional, which corresponds to the partition function in thermodynamics. When the system is in a certain 
phase, the generating functional is usually analytic for some choice of parameters, such as the current masses 
of the fermions, the chemical potential, and the temperature. The generating functional often exhibits a non- 
analytic character while the phase transitions occur. So the location and characteristics of the chiral phase 
transition in the system can be determined by the behaviors of this quantity with respect to the corresponding 
parameters. In this case, a phase transition, in which the first order derivative of the generating functional 
with some of the parameters is discontinuous, is referred to as first order or discontinuous phase transition. 
Second order or continuous phase transition exhibits continuity in first order derivative and a discontinuity 
or infinity in second order derivative. Nevertheless, thanks to the non-perturbative properties of QCD in 
the low energy region, such discussions are difficult and complicated. Accordingly, researchers often resort 
to some effective models, for example, see Refs. [43, 163, 164, 165, 166, 167, 168, 169, 170, 171, 172]. In 
the following, we take the NJL model, which is successful as well as popular in chiral phase transition 
studies [173, 174, 175, 176, 177, 178, 179, 161, 180, 181, 182], as an example, and give some discussions 
on the applications of some susceptibilities [36], similar discussions with in the framework of QED 3 can be 
found in Ref. [183]. For the details of the NJL model, please see Refs. [173, 184, 185, 186]. 

The Lagrangian of the NJL model is (here we still work in the Euclidean space, and take the number of 
flavors Nf = 2 while the number of colors Nc = 3) 

-Cnjl = ^0+ g^int m)q + g[{qqf + (gfqs'rg)^], (113) 

where g is a coupling constant with the dimension of mass“^, and the flavor and color indices are suppressed. 
Now we introduce the definitions of four kinds of susceptibilities: the chiral susceptibility Xs, the quark 
number susceptibility Xg, the vector-scalar susceptibility Xvs, and another auxiliary susceptibility Xm- For 
mathematical convenience, we first introduce these susceptibilities in the free quark gas case, where the 
interaction term in the Lagrangian is zero, £int = 0. Denoting them with the superscript their definitions 
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and expressions are then 


( 0 ) ^ _ 


Xs = 


^( 0 ) 

A,q 


d{qq)f _ NcNf 


y(0) 


do) 


dm 

d{qU)f 

dfj, 

d{qq)f 


NrNf 


/o 

pA 


-^s(rt + -gsm 


dp, 


p‘^l3gip)dp, 


Xm — 


d{qU) 


dm 


NcNf f mp^P 
TT^ Jo Eq 

- = 

A.VS 1 


h{p)dp, 


(114) 

(115) 

(116) 
(117) 


where A is a three-momentum cutoff which is introduced to avoid the ultraviolet divergence in the NJL 
model, P = 1/T, Eq = ^/mPT^, g{p) + h{p) = 2n(/r)(l - n(/r)), g{p) - h{p) = 2m(/r)(l - m{p)), 
f{p) = 1 - '^(m) - Mp), and 


n{m,p,p) = 
m{m,p, p) = 


1 


1 -h exp [P{Eo - p)] ’ 

1 


1 -h exp [P{Eo + p)]' 

3 ms. It should be no 

analytical expression, which is reasonable from the viewpoint of statistical mechanics. 


(118) 

(119) 


the subscript / represents the free quark gas systems. It should be noted that Xm^ and yi? have the same 


with Zf the QCD partition function in the free quark gas case. 

In the interacting case, these susceptibilities are coupled with each other. 


Xs = 

- xi°^(7‘)(l + 2gXs) ^ X^vJMXm, 

(121) 

Xq = 

-‘2gXv}{p)xvs+x'fHg-){^ j^Xq), 

(122) 

Xvs = 

-‘^gx^°Hg)xvs + x^vJig){^ j,^Xq), 

(123) 

Xm = 

dm ~ N 

(124) 


Using the iterative method, we can obtain the numerical results of these susceptibilities. For example, 
the vector-scalar susceptibility is the response of the effective quark mass (M = m — 2g{qq)) to the chemical 
potential p, and its results are shown in Fig. 7. The parameters used in this section are taken from Ref. [185], 
namely, m = 5.5 MeV, g = 5.074 GeV“^, and A = 631 MeV. 

It can be seen from Fig. 7 that, when T is smaller than a critical value Tc = 35 MeV, there always exists 
a convergent discontinuity of Xvs, corresponding to a first order phase transition; when T = Tc, Xvs displays 
a sharp and narrow divergent peak, which implies a second-order phase transition, or in other words, here 
is a CEP on the phase diagram; when T > T^, the discontinuity disappears and a rather broad peak of 
finite height is shown, corresponding to the crossover region. Comparing this with the results of the chiral 
susceptibility Xs and the quark number susceptibility Xq, shown in Fig. 8, we can conclude that at low 
temperature, the first order phase transition occurs at almost the same chemical potential; while in the 
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Figure 7: Xvs at different fi and T, taken from Ref. [36]. 


crossover region, if we pick the peak of these susceptibilities as the artificial critical points, we see that they 
tend to occur at different chemical potentials as the temperature increases. 

The calculated results of Xvs, Xs, and Xq in the crossover region are shown in Figs. 9 and 10, respectively. 
We can find that they exhibit different behaviors: the chiral susceptibility Xe exhibits an obvious band, so it 
is convincing to define the peak of Xs as the artificial critical points; in the high T and/or low /r region, the 
vector-scalar susceptibility Xvs tends to vanish; while the global shape of the quark number susceptibility Xq 
is just similar to the ones of Xs and Xvs, but it is non-vanishing in the high T and/or high /i region whose 
behavior is closely linked to the quark number density. Therefore, Xq can not describe the crossover property 
well in the high T and/or high /i region. 

It is very interesting and meaningful to compare the results of these susceptibilities with that of the 
thermal susceptibility xt = d{qq)/dT. For mathematical convenience, we also define Xn = d{q^q)/dT. By 
the same process employed above, we obtain a set of coupled equations for xt and Xn as 

,d{qq). 


_ d{qq) d{qq) dM d{qq) 
XT = + 


dT 


dM dT 


V )l 


dfi dT ^ dT 


=2ffxi°H/^)XT - ^xi?(M)Xn + Mf3x^°\n) - ^J-f3Xmi^J-)y 


_ d{qh) _ d{qU) dM d{qU) dfi 

arr Q 71 ,r ' a.. arr ' V ) 


dT 


dM dT dfj. dT 


=‘^9XmMxT - -^X^qXiAXn - M/3X9°Hm) 


dT ^ 

N,Nf 


(125) 


(126) 


p^El3'^h{n)dp, 


with i?= M^ + . The behavior of xt in the crossover region is shown in Fig. 10, which is very similar to 

that of Xs- 

Now we compare the results of different susceptibilities in a chiral phase diagram of QCD, as presented in 
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Fig. 11. The corresponding lines in the crossover region are determined by the peaks Xvs-, XT, and Xs- Here 
it should be stressed that, actually there is a very large number of studies investigating the phase diagram 
together with the CEP of QCD, and in various models [52, 161, 162, 171, 187, 188, 189, 190, 191, 192, 193] as 
well as functional and Lattice QCD studies [23, 194, 195, 196, 197, 198, 199, 200] the results vary quite widely, 
most of them Hnd a higher value for Tc and a lower value for fj,c of the CEP we plotted here, as summarised in 
Fig. 12 (here /is = 3/i), which is taken from Fig. 4 of Ref. [52]. The points in Fig. 12 represent effective model 
and Lattice QCD predictions, the two dashed lines are parabolas with slopes corresponding to Lattice QCD 
calculations, and the red circles are locations of the freeze-out points for heavy ion collisions at corresponding 
center of mass energies per nucleon (indicated by labels in GeV). It is shown clearly that, different models 
or parameter sets give quite different predictions (please refer to Ref. [52] for more details of these points 
and lines). The parameters we adopted here (NJL89a in Fig. 12) are just for the convenience of showing the 
qualitative results of different susceptibilities more explicitly. For a recent parameter-independent attempt 
of such studies, please see Ref. [201]. 

We see from Fig. 11 that these susceptibilities split when the system tends to lower quark chemical 
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Figure 10: Xq XT in the crossover region, taken from Ref. [36]. 


potential and/or higher temperature. As we discussed above, up to now there is no exact order parameter 
for studies of QCD chiral phase transition when beyond the chiral limit, and different susceptibilities would 
characterize different physical properties of the system. Hence, related topics deserve further discussions. 


5. Summary 

In the Dyson-Schwinger equations (DSEs) framework of QCD, we summarize studies of various vacuum 
susceptibilities together with some of their applications to the chiral phase transition of QCD. In Section 1, 
we give a brief introduction of the vacuum susceptibilities together with their related applications. Then a 
general derivation of the vacuum susceptibility using the QCD sum rules external field formula is explicated 
in the following Section 2, as well as a model independent expression, which is expressed with the quark 
propagator (two-point Green function) and the corresponding vertex function (three-point Green function). 
In the next Section 3, we review the calculations of the quark propagator and the vertex function under 
the framework of DSEs-BSEs hrst, and then use these theoretical methods and results to discuss the calcu¬ 
lations of the tensor, the vector, the axis-vector, the scalar, and the pseudo-scalar vacuum susceptibilities 
further. For the tensor vacuum susceptibility, we introduce the calculation results of the Maris-Tandy gluon 
propagator model using the rainbow-ladder approximation of DSEs. For the vector and axial-vector vac¬ 
uum susceptibilities, we present the model independent results obtained with the help of the vector and 
axial-vector Ward-Takahashi identities: the vector vacuum susceptibility is strictly 0, while the axial-vector 
vacuum susceptibility equals 3/4 of the square of the pion decay constant. For the scalar vacuum suscepti¬ 
bility, we show the results using a Gauss model (infrared part of the Maris-Tandy model) gluon propagator 
with the Ball-Chiu vertex ansatz, and comparisons with those from the bare vertex (the rainbow ladder trun¬ 
cation) are also given. We also show a straightforward explanation of a mismatch between extant estimates 
of the pseudo-scalar vacuum susceptibility. These works promoted our understandings of the five kinds of 
vacuum susceptibilities. In a word, by this time only the scalar and the tenser vacuum susceptibilities are 
still mo del-dependent, and to choose a more reliable model framework then becomes very important. With 
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the help of an effective model of QCD then, in Section 4 we show the applications of some susceptibilities to 
the QCD chiral phase transition as well as some discussions of the results. 

Here it also needs to be pointed out that, in the QCD sum rules external field approach, all kinds of vac¬ 
uum condensates and vacuum susceptibilities are introduced as independent phenomenological parameters, 
and their values are determined by fitting the theoretical predictions from QCD sum rules to the experi¬ 
mental results. The authors of Ref. [202] proposed for the first time that the four-quark condensate may 
relate to the corresponding vacuum susceptibility. If this conclusion is correct, the QCD sum rules external 
field formula can then introduce less independent parameters, which will certainly increase its predictiveness. 
Last but not least, it should be noted that up to now people still do not know how to rigorously define some 
of the vacuum condensates (including but not limited to the four-quark condensate) from the first principle 
of QCD. There is no doubt that, the issues related to the four-quark condensate and the relations between 
the four-quark condensate and the corresponding vacuum susceptibility are worth further studying. And of 
course, so are many discussions in this paper. 
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